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ABSTRACT 

The most powerful blazars are the flat spectrum radio quasars whose emission is dominated by a Compton com¬ 
ponent peaking between a few hundred keV and a few hundred MeV. We selected two bright blazars, PKS 2149- 
306 at redshift z=2.345 and S5 0836+710 at z=2.172, in order to observe them in the hard X-ray band with the 
NuSTAR satellite. In this band the Compton component is rapidly rising almost up to the peak of the emission. 
Simultaneous soft-X-rays and UV-optical observations were performed with the Swift satellite, while near-infrared 
(NIR) data were obtained with the REM telescope. To study their variability, we repeated these observations for 
both sources on a timescale of a few months. While no fast variability was detected during a single observation, 
both sources were found to be variable in the X-ray band, up to 50%, between the two observations, with larger 
variability at higher energies. No variability was detected in the optical/NIR band. These data together with 
Fermi-LAH, WISE and other literature data are then used to study the overall spectral energy distributions (SEDs) 
of these blazars. Although the jet non-thermal emission dominates the SED, it leaves the UV band unhidden, 
allowing us to detect the thermal emission of the disc and to estimate the mass of the black hole. The non-thermal 
emission is well reproduced by a one-zone leptonic model. The non-thermal radiative processes are synchrotron, 
self-Compton and external Compton using seed photons from both the broad-line region (BLR) and the torus. We 
find that our data are better reproduced if we assume that the location of the dissipation region of the jet, i?diss, is 
in-between the torus, (at i?torus), and the BLR (i?torus > Rdiss > Rblr)- The observed variability is explained 
by changing a minimum number of model parameters by a very small amount. 

Subject headings: galaxies: active - quasars: general - X-rays: general — blazar: individual: PKS 2149-306, S5 
0836+710 


1. INTRODUCTION 

Blazars are a sub-class of AGN whose emission is domi¬ 
nated by a relativistic jet pointing toward us. Their emission 
extends from the radio band to the 7 -ray and TeV band, and 
is dominated by non-thermal processes. Their spectral energy 
distribution (SED) is dominated by two humps. The first, span¬ 
ning the infrared to the X-ray band, is usually attributed to syn¬ 
chrotron emission, while the second, located from the X-ray to 
the 7 -ray band, is attributed to the inverse Compton scattering 
process. This second component can be produced by energetic 
electrons scattering only their own synchrotron photons (Syn¬ 
chrotron Self-Compton, SSC for short), or also scattering ra¬ 
diation produced externally to the jet (i.e. External Compton, 
EC for short). The latter process is likely to be important es¬ 
pecially in those sources for which the second component is 
largely dominating the SED. This usually occurs for flat spec¬ 


trum radio quasars (ESRQ), which are among the most lumi¬ 
nous persistent sources of the Universe. As such, they can be 
detected in almost all bands also at high redshifts, well above 
z > S (Romani et al. 2004, Sbarrato et al. 2012). 

The most powerful ESRQ are therefore very bright hard X- 
ray and 7 -ray sources, with the hard X-ray band (20-150 keV) 
more effective in selecting bright ESRQ at z > 4 (Ajello et al. 
2009, Ghisellini et al. 2010a), while the 7 -ray band (0.1-10 
GeV) is very effective up to z = 2 — 3 (Ghisellini et al. 2011, 
Shaw et al. 2013, Giommi et al. 2013). Hard X-ray and 7 -ray 
surveys recently obtained with the Burst Alert Telescope (BAT, 
onboard the Swift satellite, Gehrels et al. 2004) and the Large 
Area Telescope (LAT, on board Fermi, Atwood et al. 2009) in¬ 
struments (Ajello et al. 2009, Cusumano et al. 2010, Nolan et 
al. 2012, Ajello et al. 2014) provide us the opportunity to see 
if the jet properties evolve with cosmic time (Volonteri et al. 
2011; Ghisellini et al. 2011, 2013). Moreover, the non-thermal 
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SED of these powerful sources leaves the UV band unhidden 
(the synchrotron hump peaks at much smaller frequencies, and 
the high-energy hump at much larger energies; Sbarrato et al. 
2013), and at these frequencies the thermal accretion flux can 
emerge and be detected. This allows us to study the relation 
between the jet power and the accretion luminosity. 

It is clear from the above that the hard X-ray and 7 -ray bands 
are crucial to characterize the properties of these sources. In 
particular, the hard X-ray band (10-100 keV) is where the 
Compton component is rapidly rising, but it is also the band 
in which up to now we did not have very sensitive instruments 
(e.g. compared to those at softer X-rays). Thanks to the advent 
of the Nuclear Spectroscopic Telescope Array (NuSTAR) satel¬ 
lite (Harrison et al. 2013), it is now possible to obtain detailed 
X-ray spectra in the 3-79 keV band for moderately bright X-ray 
sources. Therefore we selected two bright blazars, PKS 2149- 
306 at redshift z=2.345 and S5 0836H-710 at z=2.172, and or¬ 
ganised simultaneous observing campaign with the NuSTAR 
and Swift satellites plus on-ground optical/NIR observations for 
PKS 2149-306. Both are detected in the 7 -ray band by Fermi- 
LAT. 

Both sources are high-redshift, bright blazars with a SED 
dominated by the Compton component. PKS 2149-306 was 
detected up to 100 keV with the PDS detector on board the 
BeppoSAX satellite with a very flat hard X-ray spectrum (en¬ 
ergy spectral energy index a ~ 0.4, Elvis et al. 2000). A 
flat hard X-ray spectrum up to 100 keV was later also detected 
with INTEGRAL (Bianchin et al. 2009) and SwiftlBAT (Sam- 
bruna et al. 2007). This source is detected in the 7 -ray band by 
the Term/-LAT detector (Nolan et al. 2012). S5 0836H-710 has 
been detected up to 100 keV with the PDS detector on board the 
BeppoSAX satellite with a very flat hard X-ray spectrum (en¬ 
ergy index a ~ 0.4, Tavecchio et al. 2000). The source has also 
been detected both by INTEGRAL (Beckmann et al. 2006) and 
5w//f/BAT (Sambruna et al. 2007). In particular, the BAT de¬ 
tection shows a steeper spectrum, a ~ 0 . 8 , and a flux a factor 
of ~ 5 weaker than the one recorded with BeppoSAX (Sam¬ 
bruna et al. 2007). This blazar is a bright and variable 7 -ray 
source already detected by the Energetic Gamma-ray Experi¬ 
ment Telescope (EGRET, onboard the Compton Gamma Ray 
Observatory, Hartman et al. 1999), with flare activity recently 
seen with the Termi-LAT detector (Akyuz et al. 2013). 

Here we present the results of simultaneous observations ob¬ 
tained in the X-ray band with NuSTAR and with the X-ray 
Telescope (XRT, Burrows et al. 2005) on board Swift for the 
two blazars. Eor both sources we repeated these observations 
over a time scale of a few months to check for source vari¬ 
ability. The X-ray data are complemented with simultaneous 
optical/UV data taken with the Ultra-Violet Optical Telescope 
(UVOT, Roming et al. 2005) on board Swift and, for PKS 2149- 
306, with optical/NIR data from the robotic Rapid Eye Mount 
(REM) telescope (Zerbi et al. 2004) located in La Silla, Chile. 
Eor both sources we also analyze the Termi-LAT data over a 1 
year time scale centered on om NuSTAR and Swift observations. 
While we were in the process of submitting this work, a pa¬ 
per appeared on astro-ph presenting these data of S5 0836H-710 
(Paliya 2015). The paper of Paliya is more concentrated on the 
7 -ray variability, but it also analyses and discusses the SED of 
this source. We will briefly compare our results with the one of 
that paper in the discussion. 

In this work, we adopt a flat cosmology with Hq = 70 km 

'^Data retrieved from the WISE All-Sky Source Catalog: 


s ^ Mpc ^ and Hm = 0.3. 

2. OBSERVATIONS AND DATA ANALYSIS 

The NuSTAR, Swift (and for PKS 2149-306 also REM) ob¬ 
servations are integrated with data from the Wide-field Infrared 
Survey Explorer ( 14754 ^ satellite (Wright et al. 2010) and with 
archival data from NASA/IPAC Extragalactic Database (NED) 
and the ASI Science Data Center (ASDC). 

2.1. NuSTAR observations 

The NuSTAR satellite carries two co-aligned hard X-ray tele¬ 
scopes, each consisting of a mirror module focusing high- 
energy X-ray photons in the band 3 — 79 keV onto two inde¬ 
pendent shielded focal plane modules (EPMs), referred to here 
as EPMA and EPMB (Harrison et al. 2013). 

The NuSTAR satellite observed PKS 2149-306 on 2013 De¬ 
cember 17 (obsID 60001099002) and on 2014 April 18 (obsID 
60001099004). The total net exposure times were 38.5 ks and 
44.1 ks, respectively. S5 0836H-710 was observed by NuSTAR 
on 2013 December 15 (obsID 60002045002) and on 2014 Jan¬ 
uary 18 (obsID 60002045004), for total net exposure times of 
29.7 ks and 36.4 ks, respectively. 

The EPMA and EPMB data sets were first processed with 
the NuSTARDAS software package (v. 1.4.1) jointly developed 
by the ASI Science Data Center (ASDC, Italy) and the Cali¬ 
fornia Institute of Technology (Caltech, USA). Event files were 
calibrated and cleaned with standard Altering criteria with the 
nupipeline task using the version 20140414 of the NuSTAR 
CALDB. 

Eor all four observations the EPMA and EPMB spectra of 
the sources were extracted from the cleaned event files using a 
circle of 20 pixel (^ 49 ") radius, while the background was 
extracted from two distinct nearby circular regions of 50 pixel 
radius. The ancillary response files were generated with the 
numkarf task, applying corrections for the point spread func¬ 
tion (PSP) losses, exposure maps and vignetting. All spectra 
were binned to ensure a minimum of 20 counts per bin. Both 
sources are quite bright and well detected by NuSTAR up to 79 
keV. 

2.2. Swift observations 

Swift (Gehrels et al. 2004) observed PKS 2149-306 on 2013 
December 16-17 (obsIDs 00031404013 and 00031404014) and 
on 2014 April 18 (obsID 00031404014) and S5 0836H-710 was 
observed on 2013 December 16 (obsIDs 00080399001) and on 
2014 January 18 (obsID 00080399002). 

2.2.1. XRT observations 

The XRT on board Swift is sensitive to the 0.3 — 10 keV X- 
ray energy band (Burrows et al. 2004). All XRT observations 
were carried out using the most sensitive Photon Counting (PC) 
readout mode. The XRT data sets were first processed with the 
XRTDAS software package (v.3.0.0) developed at the ASI Sci¬ 
ence Data Center (ASDC) and distributed by HEASARC within 
the HEASoft package (v. 6.16). Event files were calibrated and 
cleaned with standard Altering criteria with the xrtpipeline task 
using the calibration files available in the version 20140709 of 
the Swift-XRT CALDB. 

The two individual XRT event files for the 2013 December 
observations of PKS 2149-306 were merged together using the 
XSELECT package for a total net exposure time of 8.0 ks. The 
net exposure time of the 2014 April observation was 6.4 ks. 


http://irsa.ipac.caltech.edu/. 
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For the two observations of S5 0836+710 the total net exposure 
times were 1.9 ks and 4.7 ks, respectively. For all observations 
the energy spectra were extracted from the summed cleaned 
event hies. Events for the spectral analysis were selected within 
a circle of 20 pixel (~ 47 ") radius, which encloses about 90% 
of the PSF, centered on the source position. The background 
was extracted from a nearby circular region of 50 pixel radius. 
The ancillary response files were generated with the xrtmkarf 
task, applying corrections for the PSF losses and CCD defects 
using the cumulative exposure map. The source spectra were 
binned to ensure a minimum of 20 counts per bin. 

2.2.2. UVOT observations 

UVOT observations were performed with all six optical and 
UV lenticular biters (namely W2, M2, Wl, u, b, v) (Roming 
et al. 2005). We performed aperture photometry for all bi¬ 
ters in all observations using the standard UVOT software dis¬ 
tributed within the HEAsoft package (version 6.15.1) and the 
calibration included in the latest release of the CALDB. Counts 
were exbacted from aperture of 5 arcsec radius for all biters 
and converted to buxes using the standard zero points (Poole et 
al. 2008). The buxes were then de-reddened using the appro¬ 
priate values of E{B — V) taken from Schlegel et al. (1998) 
and Schlaby et al. (2011) with A\/E{B — V) ratios calculated 
for UVOT biters using the mean Galactic interstellar extinction 
curve from Eitzpatrick (1999). These buxes were then used to 
build the SEDs of both sources (see Eig. |^. No variability was 
detected within single exposures in any biter. The processing 
results were carefully veribed checking for possible contamina¬ 
tion from nearby objects within the source apertures and from 
objects falling within background apertures. Both sources are 
well detected in all biters. The values of the magnitudes in the 
Vega system are given in Table [T] 

2.3. REM observations 

We observed PKS 2149-306 with the Rapid Eye Mount Tele¬ 
scope (REM, Zerbi et al. 2004), a robotic telescope located 
at La Silla Observatory (Chile). It performed photometric ob¬ 
servations in the optical gri and near-infrared (NIR) JHK bi¬ 
ters in the nights of 2013 December 15 and 19. REM has a 
Ritchey-Chretien conbguration with a 60 cm f/2.2 primary and 
an overall f/8 focal ratio in a fast-moving alt-azimuth mount that 
provides two stable Nasmyth focal stations. The two cameras, 
REMIR (Conconi et al. 2004) for the NIR and ROS2 (Molinari 
et al. 2014) for the optical, both have the same held of view 
of 10 X 10 '. The telescope is able to operate in a fully au¬ 
tonomous way (Covino et al. 2004) and data are reduced and 
analyzed following standard procedures. Aperture photometry 
was derived by means of custom tool^ and calibration was 
based on objects in the held of view reported in the APASSf^ 
and 2MASSp]catalogues in the optical and NIR, respectively. 

Table reports the observed gri (AB) and JHK (Vega) 
photometry measured on the two nights, not corrected for the 
Galactic extinction of E{B — V) = 0.02 from Schlegel et al. 
(1998). 


2.4. Fermi-LAT observations 

Both sources are also bright 7 -ray emitters and are regularly 
detected by Fermi-LAT. We analyzed the data collected for one 
year between 2013 June 01 - 2014 June 01 (MJD 56444 - 
56809) following the standard procedur^^ using the Fermi- 
LAT analysis software ScienceTools v9r34pl with the 
P7REP_SOURCE_V15 instrument response functions. Events 
in the energy range 100 MeV-300 GeV were extracted within 
a 15° acceptance cone of the Region of Interest (ROI) centered 
on the location of each source. Gamma-ray buxes and spec¬ 
tra were determined by an unbinned maximum likelihood bt 
with gtlike. The background model included all known 7 - 
ray sources within the ROI from the 2nd Fermi-LAT catalog 
(Nolan et al. 2012)|^ Additionally, the model included the 
isotropic and Galactic diffuse emission components. Elux nor¬ 
malization for the diffuse and background sources were left free 
in the btting procedure. 

The Fermi-LAT light-curves above 100 MeV binned over a 
time scale of seven days are shown in Eig. [T]for both sources. 
We mark with vertical dotted lines the dates of the NuSTAR ob¬ 
servations. Eor bins with test statistic TS < 4, 95% conbdence 
level upper limits are plotted (for the meaning of TS see Mattox 
et al. 1996). We btted the 1-year Fermi-LAT spectrum for each 
source with a power-law model, bnding a best-bt photon spec¬ 
tral index of F = 2.89 ± 0.09 (TS=210) and F = 2.69 ± 0.10 
(TS=284) for PKS 2149-306 and S5 0836+710, respectively 
(errors are at 68 % conbdence interval for one parameter of in¬ 
terest). The corresponding average 1-year buxes are E{> 100 
MeV)= (8.4 ± 0.9) x 10"® ph cm-^ s"! and F(> 100 
MeV)= (5.6 ± 0.8) x 10"® ph cm-^ s"!. 

In the case of S5 0836+710 we have enough statistics to ex¬ 
tract a 1-month (2014 January 1 - 2014 February 1) source 
spectrum centered around the second NuSTAR observation. The 
btting result yielded F = 2.7±0.2 (TS=46) with an average bux 
above 100 MeV of (8.7±1.9) x 10“® ph cm“^ s“^. From Fig. 
[T] it is also apparent that a couple of months later the source 
became much more active in the Fermi-LAT band. 

2.5. X-ray spectral analysis 

Eor all observations of both sources simultaneous bts of the 
XRT and NuSTAR spectra were performed using the XSPEC 
package. In both cases a broken power-law model with an 
absorption hydrogen-equivalent column density bxed to the 
Galactic value was adopted (A^h = 1.6 x 10^° cm“^ and Wh = 
2.8 X 10^° cm-2 for PKS 2149-306 and S5 0836+710, respec¬ 
tively, Kalberla et al. 2005) . To allow for the cross-calibration 
uncertainties between the three telescopes (two NuSTAR and 
one Swift), a multiplicative constant factor has been included, 
kept equal to 1 for the EPMA specba and free to vary for the 
EPMB and XRT spectra. In the case of EPMB the difference 
is in the range of 2^%, while for XRT it is somewhat larger 
but alway less than 10%. This is consistent with the values usu¬ 
ally found for other sources. We bnd that this model provides a 
good description of the observed spectra in the 0.3-79 keV en¬ 
ergy band. The results of the spectral bts are in Table while 
Eigurej^ shows an example of the XRT and NuSTAR spectra to- 


*^https://pypi.python.org/pypi/SRPAstro.FITS/ 

*®http://www.aavso.org/apass 
'^http://www.ipac.caltech.edu/2mass/ 

**http://fenni.gsfc.nasa.gov/ssc/data/analysis/ 

In order to test the effect of possible new gamma-ray sources not included in 2FGL, we reanalyzed the data using the Fermi-LAT 4-year point source catalog 
(3FGL, Acero et al. 2015), finding fully consistent results. 
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PKS 2149-306 


Date 

V 

b 

u 

Wl 

M2 

W2 

17 Dec 2013 

18 Apr 2014 

17.83 ±0.11 
17.61 ±0.05 

17.86 ±0.06 
17.75 ±0.02 

17.14 ±0.06 
17.14 ±0.02 

18.22 ±0.10 
18.06 ±0.04 

20.46 ±0.07 
20.00 ±0.11 

20.12 ±0.08 

S5 0836+710 

15 Dec 2013 
18 Jan 2014 

16.95 ±0.03 

17.12 ±0.02 

16.22 ±0.02 

16.92 ±0.03 

17.45 ±0.03 
17.53 ±0.04 

17.91 ±0.03 


Table 1 

UVOT Vega v, h, u, Wl, M2, W2 observed magnitudes oe PKS 2149-306 and S5 0836+710 (magnitudes not corrected 

EOR Galactic extinction). 


Date 

a 

r 

i 

J 

H 

Ks 

15 Dec 2013 
19 Dec 2013 

17.82 ±0.05 
17.75 ±0.05 

17.58 ±0.06 
17.36 ±0.06 

17.40 ±0.06 
17.38 ±0.06 

16.56 ±0.06 
16.58 ±0.07 

16.00 ±0.08 
15.89 ±0.07 

15.10 ±0.11 
15.22 ±0.13 


Table 2 

REM AB (gri) AND VEGA (JHK) OBSERVED MAGNITUDES OE PKS 2149-306 (magnitudes NOT CORRECTED FOR GALACTIC 

EXTINCTION). 
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Fig. 1.— Gamma-ray light curves (> 100 MeV) of S5 0836+710 (top panel) and PKS 2149-306 (bottom panel) as observed by FermifLAT. The vertical dotted 
lines mark the dates of the NuSTAR observations. Upper limits correspond to TS < 4 (see text). 
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PKS2149-306 


Date 

Ti 

r2 

break 

keV 

.F2-10kev 
erg cm“^ 

-^10—40/,kev 

ergcm“^ 

/ dof 

2013 Dec 17 

2014 Apr 18 

^•^^-0.08 

Q q7+0-07 
-0.09 

1 ok+0.02 

-L-O0_o.oi 

-^•^'^-0.02 

^■^-^-0.46 
q or-|-1.36 
*^•^^-0.63 

2.2 X 10"“ 

1.9 X 10"“ 

5.0 X 10"“ 

3.7 X 10-“ 

1095.8/ 1108 

936.2/952 

S5 0836H-710 

Date 

Ti 


break 

keV 

.F2-10kev 
erg cm“^ 

^10—40 kev 

erg cm“^ s“^ 

X^ / dof 

2013 Dec 15 

2014 Jan 18 

-L-GO_o.32 

1 1 O+0.08 

-’^•^O-O-IO 

-L-OU_o.o2 

-L-OO 0.01 

1 '7q + -L.27 

‘'^-0.48 

2-841^62 

1.6 X 10"“ 

2.5 X 10"^^ 

2.3 X 10"“ 

3.6 X 10"“ 

642.6/611 

896.4 / 892 


Table 3 

Parameters of the X-ray spectral analysis for the simultaneous fit of the NuSTAR and Swift/XRT data. The errors are 
AT 90% level of confidence for one parameter of interest. Fluxes are corrected for the galactic absorption. 


data and folded model 



Energy (keV) 

Fig. 2 .— X-ray spectrum of PKS 2149-306 as observed by 5vyr/r/XRT and NuSTAR on 2014 April 18, together with the broken power-law best-fit model. 
NuSTAR data are filled circles and triangles (black and red, in the electronic version), while SwiftfXRT data are empty squares (green in the electronic version). 
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gether with the best-fit model for the observation of PKS 2149- 
306 on 2014 April 18. To test the robustness of the spectral 
curvature, we estimated the instrumental cross-calibration fac¬ 
tors by fitting the data in a common energy band (3-9 keV) 
adopting a power-law model. We found fully consistent best- 
fit values of the cross-calibration factors. We then verified that 
the broken power-law best-fit spectral results were unchanged 
using these values. Besides using a broken power-law model 
we tried a simpler power-law model, but this did not provide a 
good fit to the data for either source. In fact, in all cases the 
reduced Xr was larger than 1.1-1.2 with more than 900 d.o.f.. 
The F-test shows that the improvement of these values obtained 
with the broken power-law model is significant with a probabil¬ 
ity smaller than 10 “^° in three cases and smaller than 10 “^^ in 
one case. Moreover, when using the power-law model the con¬ 
stant value included to take into account the cross-calibration 
uncertainties was too large or too low (e.g. 20-30% difference) 
with respect to the values usually found, another indication that 
the single power-law model is not correct. 

In Table[^we also provide the fluxes in the two bands 2 — 10 
and 10 — 40 keV for the four observations. Both sources have 
varied between the two observations, spanning a time scale of 
one month for S5 0836H-710 and of four months for PKS 2149- 
306. No fast variability is observed during the single observa¬ 
tions. However, while for S5 0836H-710 there is an increase by 
a factor of ^ 1.5 over the full 0.3-79 keV band, with a hint of 
lower variability below ~ 1 keV, for PKS 2149-306 a variabil¬ 
ity of ~ 30% is present essentially only above 10 keV. To better 
show this we plot all XRT and NuSTAR spectra in Figure 

3. OVERALL SPECTRAL ENERGY DISTRIBUTION 

Fig. |4 shows the X-ray data of both sources, and compares 
them with archival observations. For S5 0836H-710 the two ob¬ 
servations correspond to different flux levels, and slightly dif¬ 
ferent spectra. The XRT spectrum smoothly joins the NuSTAR 
data points in both observations. The two data sets differ mainly 
around 10 keV, by a factor ~1.6, with less variations at both 
ends of the observed X-ray band. On the contrary, for PKS 
2149-306 the variability amplitude monotonically increases at 
higher frequencies, while the two spectra are very similar in the 
soft X-ray band. The archival data (orange symbols) indicate 
that both sources (but especially S5 0836H-710) can vary in X- 
rays by about an order of magnitude. Note also the much softer 
archival spectrum of PKS 2149-306 by XMM (and by INTE¬ 
GRAL). 

Fig. 0 shows the overall SED of the two sources. Both of 
them show, besides the typical double-peak SED of blazars, a 
third narrow peak in the IR-UV band, which we identify as 
thermal emission from the accretion disk. Eitting it with a stan¬ 
dard Shakura & Sunyaev (1973) disk we find both the disk lu¬ 
minosity Ld and the black hole mass M (as listed in Table 
and Table see Calderone et al. 2013 for a full discussion 
about this method). Eor a given efficiency p (defined through 
Ld = r]Mc^), and when the peak of the disk emission is visi¬ 
ble, this method returns a value of the black hole mass with a 
relatively small uncertainty (factor ~1.5), better than the virial 
method (factor 3^, Vestergaard et al. 2006; Park et al. 2012). 
Adopting rj = 0.08, we find a black hole mass of 5 x 10 ®Mq for 
S5 0836H-710) and 3.5 x 10®M© for PKS 2149-306j3 These 
values differ slightly from what we reported previously using 


the same method (in Ghisellini et al. 2010a, hereafter GGIO). 
This is due to the better coverage of the IR part of the spec¬ 
trum now available from WISE and REM. With these values of 
the black hole mass, the disks of both sources emit at ~ 20 - 
30% of the Eddington luminosity. The blazar S5 0836H-710 has 
also been monitored for reverberation mapping by Kaspi et al. 
(2007). Using the observed lag between the continuum, the 
CIV EWHM and Eq. 5 of Kaspi et al. (2000), they obtained a 
black hole mass M = 2.6 x 10®M©. Instead, using the UV lu¬ 
minosity, the CIV EWHM and Eq. 7 of Vestergaard & Peterson 
(2006), they derive M = 1.8 x 10^°Mq. 

As discussed also in GGIO, these sources are weak in the 7 - 
ray band. This is due to two effects. The first is a k-correction 
effect; for increasing redshifts, the observed hard X-ray spec¬ 
trum is closer to the high-energy peak, and thus brighter. The 
second effect is due to the intrinsic shift of the high-energy 
peak frequency as the bolometric non-thermal luminosity is in¬ 
creased. This implies that the most powerful jetted sources can 
be better found in hard X-ray surveys rather than in 7 -ray sur¬ 
veys (GGIO). Eor S5 0836 h- 7 10 Fermi-LAT detected the source 
in one-month integration times around the second NuSTAR ob¬ 
servation (filled red points with black circles) while PKS 2149- 
306 is detected using one-year integration times (from June 1st 
2013 to June 1st 2014). Both 7 -ray spectra are steep, as usual 
for powerful ESRQ (see e.g. Nolan et al. 2012). The fact that 
the synchrotron far-mm-IR spectrum is steep is consistent with 
the assumption that both the low and high-energy non-thermal 
peaks of the SED are produced by the same population of elec¬ 
trons. 

We applied a one-zone leptonic model, fully described in 
Ghisellini & Tavecchio (2009), to interpret the SED of the two 
sources, and to explore the variability that both sources experi¬ 
enced between the two Swift-i-NuSTAR observations. In brief, 
the model assumes that the emitting source is a homogeneous 
sphere located at a distance i?diss from the black hole, in a coni¬ 
cal jet of semi-aperture angle (j) = 0.1 rad. The BLR is assumed 
to be a spherical shell of radius TJblr = cm and the 

size of the IR torus is assumed to be i?torus = 2 x 
cm (here Ld ,45 is the disk luminosity in units of 10 ^® erg s“^). 
The source is assumed to move with a bulk Lorentz factor F 
in a direction making an angle 6^ with the line of sight. The 
magnetic field B (as measured in the comoving frame) depends 
upon the distance from the black hole and the bulk Lorentz fac¬ 
tor, to yield a constant Poynting flux Pb oc 

The energy distribution of the emitting particles is found 
through the continuity equation, accounting for continuous in¬ 
jection, radiative cooling, and electron-positron pair produc¬ 
tion. Electrons (with random Lorentz factors 1 < 7 < 7 max) 
are assumed to be injected with a total power P[ (in the co¬ 
moving frame) throughout the source with a broken power-law 
distribution, cx 7 “'*! and cx 7 “®= below and above 7 b, respec¬ 
tively. The radiative processes are synchrotron, self-Compton 
and external Compton with both photons from the BLR and the 
torus. We also consider the presence of an X-ray corona, emit¬ 
ting 30% of Ld with a spectrum F^, cx exp(—/iu/150 keV). 

Due to the compactness of the source, required to account for 
the fast variability, the synchrotron emission is self-absorbed 
at radio frequencies, up to hundreds of GHz (observer frame). 
Therefore these models cannot reproduce the radio data. 


^®The values of the black hole masses are somewhat dependent upon rj: the larger 7 , the greater the black hole mass. For instance, adopting rj = 0.06 (0.15), we 
would obtain M = 3 X lO^Mg (6 X 10®Mq) for S5 0836+710. 
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PKS 2149-306 
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Fig. 3.— X-ray SED of S5 0836+710 (left panel) and of PKS 2149-306 (right panel) as observed by Swift/XKT and NuSTAR on four different dates (see text). 
XRT data points are filled-empty circles (green and purple in the electronic version), NuSTAR/FPMA data points are filled-empty squares (red and blue in the 
electronic version), M(S7]4R/FPMB data points are filled-empty triangles (black and magenta in the electronic version). 


Name 

[1] 

M 

[2] 

r 

[3] 

0. 

[4] 

-^diss 

[5] 

R 

[6] 

-Rblr 

[7] 

P( 

[8] 

B 

[9] 

7b 

[10] 

Tmax 

[11] 

Si 

[12] 

S2 

[13] 

7c 

0836+710H 

5e9 

16 

3 

1950 

195 

1500 

0.11 

1.11 

250 

5e3 

1.7 

3.2 

5.7 

0836+710L 

5e9 

16 

3 

2100 

210 

1500 

0.09 

1.04 

190 

4e3 

1.6 

3.2 

9.8 

0836+710 GGIO 

3e9 

14 

3 

540 

54 

1500 

0.22 

3.28 

90 

2e3 

-1 

3.6 

2.1 

2149-306 H 

3.5e9 

14 

3 

1365 

137 

1212 

0.2 

1.05 

75 

4e3 

0.5 

3.3 

2.9 

2149-306 L 

3.5e9 

14 

3 

1365 

137 

1212 

0.1 

1.05 

50 

4e3 

1 

3.0 

2.9 

2149-306 GGIO 

5e9 

15 

3 

1200 

120 

1224 

0.18 

1.12 

60 

3e3 

0 

3.3 

1 


Table 4 

Input parameters used to model the SED. “L” and “H”stand for low and high state, respectively; GGIO stands for 
Ghisellini et al. (2010a), where the two sources were also studied. Col. [1]: Source name and state/observation; 
Col. [2]: Black hole mass in solar mass units; Col. [3]: bulk Lorentz factor; Col. [4]: viewing angle (degrees); Col 
[5]: DISTANCE OF THE BLOB FROM THE BLACK HOLE IN UNITS OF 10^® CM; COL. [ 6 ]: SOURCE SIZE IN UNITS OF 10^® CM; COL. [7]: 
SIZE OF THE BROAD LINE REGION IN UNITS OF 10^® CM; COL. [ 8 ]: POWER INJECTED IN THE BLOB CALCULATED IN THE COMOVING 
FRAME, IN UNITS OF lO'^® ERG S“^; COL. [9]: MAGNETIC FIELD IN GAUSS; COL. [10] AND [11]: MINIMUM AND MAXIMUM RANDOM 
Lorentz factors of the injected electrons; Col. [12] and [13]: slopes of the injected electron distribution [< 3 ( 7 )] 
BELOW and above 7 b; Col. [13]: value of the minimum random Lorentz factor of the electrons cooling in R/c. The 
SPECTRAL SHAPE OF THE CORONA IS ASSUMED TO BE OC exp(— /lu/150 keV). 
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Name 

log id 

id/iEdd 

logPr 

logPe 

logPe 

logPp 

[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

0836+710H 

47.4 

0.3 

46.8 

46.7 

45.9 

48.6 

0836+710L 

47.4 

0.3 

46.7 

46.7 

45.9 

48.5 

0836+710 GGIO 

47.4 

0.5 

46.6 

46.4 

45.5 

48.0 

2149-306 H 

47.2 

0.3 

47.1 

46.2 

45.7 

48.3 

2149-306 L 

47.2 

0.3 

46.7 

46.2 

45.6 

48.3 

2149-306 GGIO 

47.2 

0.2 

46.6 

46.2 

45.3 

48.0 


Table 5 

Accretion and jet powers: Col. [1]: Source name and state/observation; Col. [2]: Logarithm of the accretion disk 
LUMINOSITY (UNITS ARE ERG S“^); COL. [3]: ACCRETION DISK LUMINOSITY IN EDDINGTON UNITS COL. [4]-[7] LOGARITHM OF THE 
JET POWER IN THE FORM OF RADIATION (P^), POYNTING FLUX (Pb), BULK MOTION OF ELECTRONS (Pe) AND [7] PROTONS (Pp 
(ASSUMING ONE COLD PROTON PER EMITTING ELECTRON). UNITS ARE ERG S“^. THESE VALUES REFER TO ONE JET ONLY. 


E [keV, obs] 

1 10 100 



E [keV, obs] 


1 10 100 



Fig. 4.— X-ray SED of S5 0836-1-710 (left panel) and PKS 2149-306 (right panel) as observed by Swift/XKT and NuSTAR, compared with previous observations 
by Swift (green dots), INTEGRAL, XMM and BeppoSAX (grey), and other ai'chival obseiwations (orange). Our XRT and NuSTAR data points for the two observing 
periods are labelled. The solid lines refer to the model used to explain the overall SED, as shown in Fig. 
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Fig. 5.— Overall SED of S5 0836+710 (left panel) and of PKS 2149-306 (right panel). The indicated REM, UVOT, XRT and NuSTAR data are simultaneous. 
For S5 0836+710, the Fermi-LAT 7-ray data (red points encircled by black circles) are time integrated over a month (2014 Jan. lst-2014 Feb. 1st), centered on the 
second NuSTAR observation, while the bow-tie represents the result of the one-year integration (2013 June lst-2014 June 1st). Fot PKS 2149-306, the Fermi~LA^ 
data (bow tie) are the result of one-year integration (2013 June lst-2014 June 1st). The solid lines refer to the one-zone leptonic model discussed in the text. The 
dashed black line is the contribution of the accretion disk, the IR torus, and the X-ray corona. We also show archival data (empty circles, orange in the electronic 
version). 


Fig. 0 shows the models corresponding to the two states of 
each source and Table |^lists the model parameters. It also re¬ 
ports the parameters used in GGIO to ht another data set. To 
model the two observed states, we changed a minimum number 
of parameters. Both sources have i?torus > Rdiss > Rblr- 
This choice is preferred (with respect to Rdiss within the BLR) 
because of the small value of the high-energy peak, requiring 
seed photons of frequency smaller than the hydrogen Lya pho¬ 
tons (that dominate the BLR emission). In this respect, NuS¬ 
TAR is crucial because the hardness of its spectrum, together 
with the extrapolation from the 7 -ray energies, greatly helps in 
pinpointing the peak frequency of the high-energy hump. This 
is the major difference with respect to the parameters reported 
by GGIO for S5 0836H-710, which assumed Rdiss < f?BLR. 

The observed variability of S5 0836H-710 is ascribed to a 
small difference in the position of the emitting region, in turn 
implying a small change in the magnetic field. Furthermore, 
the injected power changes by 10 % and the break energy of the 
injected electron distribution changes by 20 %. 

For PKS 2149-306 the variability can be explained by chang¬ 
ing the injected power by a factor of 2, and by a 50% change in 
the break energy of the injected electron distribution. 

Comparing these changes with the distribution of the param¬ 
eters in many Fermi-LAJ' blazars (as listed e.g. in Ghisellini 
et al. 2010 b), we conclude that the changes required to explain 
the observed (factor ~ 1.5) variability are very small. 

4. DISCUSSION AND CONCLUSIONS 

The simultaneous observations of NuSTAR and Swift/XRT 
revealed that the X-ray spectrum from ~0.3 to 60 keV of both 
sources are well described by a broken power-law model. Both 
indices are very hard, with no requirement of absorption in ex¬ 
cess of the Galactic one. This broken power-law behavior is 
well reproduced by the model as due to the combination of 
two effects. First, the electron energy distribution, below the 


cooling energy 7 c, retains the slope of the injected electrons 
(i.e. 7“'^0. which is very hard (see Table |^. Second, what 
we observe as low-energy X-ray emission at the frequency 
is seen, in the comoving frame, as ultraviolet emission (i.e. 

-\- z)/5), produced by the scattering between low- 
energy electrons and seed photons with frequencies smaller 
than v'^. In both sources, the main contributions to the seed 
photons come from the IR torus, the BLR, and the disk. All 
these components have a corresponding radiation energy den¬ 
sity peaking at some frequency, which is seen between the opti¬ 
cal and the UV in the comoving frame of the jet. Inverse Comp¬ 
ton scattered photons at, for example, v' ~ 10^^ Hz can be 
produced using seed photons of frequencies smaller than 10 ^^ 
Hz, which do not correspond to the full energy density of the 
seed photons. Scattered photons of higher energies, instead, can 
use the entire seed photon distribution. In other words, there is 
a paucity of seed photons at low energies, making the Inverse 
Compton spectrum harder below a few keV. What is observed in 
the two blazars studied here can then be regarded as a proof of 
the “External Compton” process, where the major contribution 
to the seed photons comes from radiation produced externally 
to the jet. Furthermore, the spectrum between 0.3 and 60 keV is 
smooth, indicating the absence of other emitting components. 

Both blazars show variability, but not extreme variability. 
This agrees with a relatively large source size, and a corre¬ 
spondingly relatively large distance from the black hole. The 
observed light crossing time predicted by our model is = 
(i?/c)(l -f z)/5, corresponding to ~13 days for S5 0836H-710 
and to ~10 days for PKS 2149-306. This agrees with no vari¬ 
ability observed within a single epoch observation, and is con¬ 
sistent with the emission site being between the BLR and the 
torus (i?BLR < 7?diss < 7?torus). We note that, contrary to what 
we have observed, significant variability on shorter time-scales 
has been reported for S5 0836H-710 by Akyuz et al. (2013). In 
the framework of our modeling this is explained by assuming 
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that the component dominating the emission at a given time is 
not always located in the same place along the jet. Sometimes 
it can be at a position closer to the black hole, and be more 
compact, giving rise to variability on shorter time-scales and to 
a slightly different SED. 

We have shown that the observed moderate variability can be 
produced by a rather small change in the injected power, break 
energy of the electron distribution and in the location of the 
emitting region. The largest change is needed for PKS 2149- 
306, requiring a factor of two change in the injected power. 

The NuSTAR data at high energies, coupled with the 7 -ray 
data, allow us to determine where the high-energy hump of the 
SED peaks. This is at ~1 MeV (observed, thus ~3 MeV, rest 
frame) for both sources. As already mentioned, this helps to 
determine the location of the emitting region, making us pre¬ 
fer an emitting region located between the BLR and the torus. 
We stress that this may not none thiecessarily always be the 
case since the emitting site can change, and sometimes it can 
be within the BLR. Eor a given electron distribution, the high- 
frequency peak would shift to higher values, in this case. The 
data of S5 0836H-710 have been analysed also by Paliya (2015), 
who fit the data using the same one single zone leptonic emis¬ 
sion model. We find consistent results, but our black hole mass 
is slightly larger, and our emitting region is just outside the 
BLR, while in Paliya et al. (2015) is inside. We find that while 
the parameters of Paliya et al. (2015) indeed satisfactorily fit 
the overall SED and the average NuSTAR spectrum, our fit ac¬ 
counts for the two states of the source in a better way, and in 
particular the two sets of SwiftfXRT+NuSTAR data. 

Both blazars have a black hole with mass exceeding 10^Mq 
and accretion disks emitting at about one-third of the Eddington 
rate. The power that the jet expends for producing the radiation 
we see (P^ in Table is of the same order as Ld (note that the 
listed values of refer to one jet only, and should be doubled). 
Pi- is a lower limit of the true jet power Pjet. An estimate of 
Pjet can be derived assuming, following Nemmen et al. (2012), 
that Pjet lOPr. Another estimate can be found assuming that 
there is one proton for each emitting electron (this is the quan¬ 
tity Pp in Table 1^. Either way, we are forced to conclude that 
the jet power exceeds the luminosity of the accretion disk, thus 
suggesting that the jet is powered not only by accretion, but also 
by the rotational energy of a spinning black hole, as found by 
Ghisellini et al. (2014) for a large sample of blazars. 

Einally, given the very hard spectrum detected by NuSTAR, 
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and their high-energy peak at ^1-10 MeV, one can wonder if 
powerful blazars can significantly contribute to the X-ray back¬ 
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